source('../settings/settings.R')
source('commonFunctions.R')
persons <- SELECTED_SUBJECTS

all_Drive2 <- read.csv('..//data/processed/distancewise/TT1_Drive_2_15m_15m.csv')
all_Drive2$Subject <- as.factor(all_Drive2$Subject)
all_Drive2$logPerspiration <- log(all_Drive2$Perspiration)

Acceleration Plot

computeThreshold <- function(x, xmin=min(x), xmax=max(x)) {
  x <- x[x >= xmin & x <= xmax]
  xArr <- matrix(x, nrow = 1, ncol = length(x))
  threshold <- otsu(xArr, range=c(xmin, xmax))
  return(threshold)
}
computeThresholdForAcceleration <- function(acc){
  return(computeThreshold(acc, xmin=0, xmax=20))
}
findModes <- function(x) {
  ux <- unique(x)
  ux[which.max(tabulate(match(x, ux)))]
}
acc_thresholds <- vector(mode="list", length=length(persons)) 
names(acc_thresholds) <- persons

acc_sd_thresholds <- vector(mode="list", length=length(persons)) 
names(acc_sd_thresholds) <- persons

speed_sd_thresholds <- vector(mode="list", length=length(persons)) 
names(speed_sd_thresholds) <- persons

for (p in persons) {
  pData <- all_Drive2[(all_Drive2$Subject==as.integer(p) | all_Drive2$Subject==p),]
  acc_th <- computeThresholdForAcceleration(pData$Acceleration)
  acc_th_sd <- computeThresholdForAcceleration(pData$Acc_std)
  speed_th_sd <- computeThresholdForAcceleration(pData$Speed_std)
  acc_thresholds[[p]] <- as.numeric(acc_th)
  acc_sd_thresholds[[p]] <- as.numeric(acc_th_sd)
  speed_sd_thresholds[[p]] <- as.numeric(speed_th_sd)
  
  acc_pl <- ggplot(pData, aes(x=Acceleration)) + geom_density() + geom_vline(aes(xintercept=acc_th), color="blue", linetype="dashed", size=1)
  acc_sd_pl <- ggplot(pData, aes(x=Acc_std)) + geom_density() + geom_vline(aes(xintercept=acc_th_sd), color="red", linetype="dashed", size=1)
  speed_sd_pl <- ggplot(pData, aes(x=Speed_std)) + geom_density() + geom_vline(aes(xintercept=speed_th_sd), color="red", linetype="dashed", size=1)
  
  plot(acc_pl)
  print(paste0(p, " - Threshold (Acc) =", acc_th, " & (Acc SD) =", acc_th_sd, " & (Speed SD) =", acc_th_sd))
  
}
[1] "01 - Threshold (Acc) =9.140625 & (Acc SD) =7.265625 & (Speed SD) =7.265625"
[1] "02 - Threshold (Acc) =10.8203125 & (Acc SD) =5.5859375 & (Speed SD) =5.5859375"
[1] "03 - Threshold (Acc) =9.3359375 & (Acc SD) =6.1328125 & (Speed SD) =6.1328125"
[1] "04 - Threshold (Acc) =9.296875 & (Acc SD) =6.328125 & (Speed SD) =6.328125"
[1] "05 - Threshold (Acc) =9.4921875 & (Acc SD) =5.4296875 & (Speed SD) =5.4296875"
[1] "06 - Threshold (Acc) =10.3515625 & (Acc SD) =6.4453125 & (Speed SD) =6.4453125"
[1] "07 - Threshold (Acc) =10.5078125 & (Acc SD) =6.4453125 & (Speed SD) =6.4453125"
[1] "09 - Threshold (Acc) =10.4296875 & (Acc SD) =7.3828125 & (Speed SD) =7.3828125"
[1] "12 - Threshold (Acc) =9.8046875 & (Acc SD) =5.9765625 & (Speed SD) =5.9765625"
[1] "13 - Threshold (Acc) =10.625 & (Acc SD) =6.6796875 & (Speed SD) =6.6796875"
[1] "15 - Threshold (Acc) =9.0234375 & (Acc SD) =7.2265625 & (Speed SD) =7.2265625"
[1] "16 - Threshold (Acc) =9.921875 & (Acc SD) =7.6953125 & (Speed SD) =7.6953125"
[1] "17 - Threshold (Acc) =10.3515625 & (Acc SD) =7.1484375 & (Speed SD) =7.1484375"
[1] "18 - Threshold (Acc) =10.5859375 & (Acc SD) =8.3984375 & (Speed SD) =8.3984375"
[1] "22 - Threshold (Acc) =10.6640625 & (Acc SD) =5.5859375 & (Speed SD) =5.5859375"
[1] "24 - Threshold (Acc) =9.296875 & (Acc SD) =5.5078125 & (Speed SD) =5.5078125"
[1] "29 - Threshold (Acc) =6.796875 & (Acc SD) =5.5859375 & (Speed SD) =5.5859375"
[1] "30 - Threshold (Acc) =9.9609375 & (Acc SD) =6.71875 & (Speed SD) =6.71875"
[1] "31 - Threshold (Acc) =11.40625 & (Acc SD) =6.6015625 & (Speed SD) =6.6015625"
[1] "32 - Threshold (Acc) =7.0703125 & (Acc SD) =5.4296875 & (Speed SD) =5.4296875"
[1] "41 - Threshold (Acc) =13.203125 & (Acc SD) =5.2734375 & (Speed SD) =5.2734375"

mean_pp <- vector(mode="list", length=length(persons)) 
names(mean_pp) <- persons

std_pp <- vector(mode="list", length=length(persons)) 
names(std_pp) <- persons

# Mean (Turning)
mean_pp_seg0 <- vector(mode="list", length=length(persons)) 
names(mean_pp_seg0) <- persons
mean_pp_seg0_1 <- vector(mode="list", length=length(persons)) 
names(mean_pp_seg0_1) <- persons
mean_pp_seg0_2 <- vector(mode="list", length=length(persons)) 
names(mean_pp_seg0_2) <- persons
mean_pp_seg0_3 <- vector(mode="list", length=length(persons)) 
names(mean_pp_seg0_3) <- persons
mean_pp_seg0_4 <- vector(mode="list", length=length(persons)) 
names(mean_pp_seg0_4) <- persons

# Mean (Straight)
mean_pp_seg1 <- vector(mode="list", length=length(persons)) 
names(mean_pp_seg1) <- persons
mean_pp_seg2 <- vector(mode="list", length=length(persons)) 
names(mean_pp_seg2) <- persons
mean_pp_seg3 <- vector(mode="list", length=length(persons)) 
names(mean_pp_seg3) <- persons
mean_pp_seg4 <- vector(mode="list", length=length(persons)) 
names(mean_pp_seg4) <- persons
mean_pp_max <- vector(mode="list", length=length(persons)) 
names(mean_pp_max) <- persons

# SD (Turning)
std_pp_seg0 <- vector(mode="list", length=length(persons)) 
names(std_pp_seg0) <- persons
std_pp_seg0_1 <- vector(mode="list", length=length(persons)) 
names(std_pp_seg0_1) <- persons
std_pp_seg0_2 <- vector(mode="list", length=length(persons)) 
names(std_pp_seg0_2) <- persons
std_pp_seg0_3 <- vector(mode="list", length=length(persons)) 
names(std_pp_seg0_3) <- persons
std_pp_seg0_4 <- vector(mode="list", length=length(persons)) 
names(std_pp_seg0_4) <- persons

# SD (Straight)
std_pp_seg1 <- vector(mode="list", length=length(persons)) 
names(std_pp_seg1) <- persons
std_pp_seg2 <- vector(mode="list", length=length(persons)) 
names(std_pp_seg2) <- persons
std_pp_seg3 <- vector(mode="list", length=length(persons)) 
names(std_pp_seg3) <- persons
std_pp_seg4 <- vector(mode="list", length=length(persons)) 
names(std_pp_seg4) <- persons
std_pp_max <- vector(mode="list", length=length(persons)) 
names(std_pp_max) <- persons

# Acceleration Partitioning
mean_pp_HighAcc <- vector(mode="list", length=length(persons)) 
names(mean_pp_HighAcc) <- persons
mean_pp_LowAcc <- vector(mode="list", length=length(persons)) 
names(mean_pp_LowAcc) <- persons

std_pp_HighAcc <- vector(mode="list", length=length(persons)) 
names(std_pp_HighAcc) <- persons
std_pp_LowAcc <- vector(mode="list", length=length(persons)) 
names(std_pp_LowAcc) <- persons

# Accel + Segmentation
mean_pp_HighAcc1 <- vector(mode="list", length=length(persons)) 
names(mean_pp_HighAcc1) <- persons
mean_pp_LowAcc1 <- vector(mode="list", length=length(persons)) 
names(mean_pp_LowAcc1) <- persons
std_pp_HighAcc1 <- vector(mode="list", length=length(persons)) 
names(std_pp_HighAcc1) <- persons
std_pp_LowAcc1 <- vector(mode="list", length=length(persons)) 
names(std_pp_LowAcc1) <- persons

mean_pp_HighAcc2 <- vector(mode="list", length=length(persons)) 
names(mean_pp_HighAcc2) <- persons
mean_pp_LowAcc2 <- vector(mode="list", length=length(persons)) 
names(mean_pp_LowAcc2) <- persons
std_pp_HighAcc2 <- vector(mode="list", length=length(persons)) 
names(std_pp_HighAcc2) <- persons
std_pp_LowAcc2 <- vector(mode="list", length=length(persons)) 
names(std_pp_LowAcc2) <- persons

mean_pp_HighAcc3 <- vector(mode="list", length=length(persons)) 
names(mean_pp_HighAcc3) <- persons
mean_pp_LowAcc3 <- vector(mode="list", length=length(persons)) 
names(mean_pp_LowAcc3) <- persons
std_pp_HighAcc3 <- vector(mode="list", length=length(persons)) 
names(std_pp_HighAcc3) <- persons
std_pp_LowAcc3 <- vector(mode="list", length=length(persons)) 
names(std_pp_LowAcc3) <- persons

mean_pp_HighAcc4 <- vector(mode="list", length=length(persons)) 
names(mean_pp_HighAcc4) <- persons
mean_pp_LowAcc4 <- vector(mode="list", length=length(persons)) 
names(mean_pp_LowAcc4) <- persons
std_pp_HighAcc4 <- vector(mode="list", length=length(persons)) 
names(std_pp_HighAcc4) <- persons
std_pp_LowAcc4 <- vector(mode="list", length=length(persons)) 
names(std_pp_LowAcc4) <- persons


for(p in persons) {
  pData <- all_Drive2[(all_Drive2$Subject==as.integer(p) | all_Drive2$Subject==p),]
  pData_act2 <- pData[pData$Activity==2,]
  
  # Data Partitioning
  pData_seg0 <- pData[pData$Phase==0,]
  pData_seg0_1 <- pData[pData$Activity==0 & pData$Phase==0 & pData$Time > 100 & pData$Time < 200,]
  pData_seg0_2 <- pData[pData$Activity==0 & pData$Phase==0 & pData$Time > 200 & pData$Time < 300,]
  pData_seg0_3 <- pData[pData$Activity==0 & pData$Phase==0 & pData$Time > 300 & pData$Time < 400,]
  pData_seg0_4 <- pData[pData$Activity==0 & pData$Phase==0 & pData$Time > 400 & pData$Time < 500,]
  
  pData_seg1 <- pData[pData$Phase==1 & pData$Activity==2 & pData$Time < 110,]
  pData_seg2 <- pData[pData$Phase==2 & pData$Activity==2 & pData$Time < 250,]
  pData_seg3 <- pData[pData$Phase==3 & pData$Activity==2 & pData$Time < 350,]
  pData_seg4 <- pData[pData$Phase==4 & pData$Activity==2,]
  
  pData_HighAcc <- pData[pData$Acceleration >= acc_thresholds[[p]],]
  pData_LowAcc <- pData[pData$Acceleration < acc_thresholds[[p]],]
  # pData_HighAcc <- pData[pData$Acc_std < acc_sd_thresholds[[p]] & pData$Speed_std < speed_sd_thresholds[[p]],]
  # pData_LowAcc <- pData[pData$Acc_std >= acc_sd_thresholds[[p]] | pData$Speed_std >= speed_sd_thresholds[[p]],]
  
  pData_HighAcc1 <- pData_HighAcc[(pData_HighAcc$Distance < 350),]
  pData_LowAcc1 <- pData_LowAcc[(pData_LowAcc$Distance < 350),]
  pData_HighAcc2 <- pData_HighAcc[(pData_HighAcc$Distance >= 350 & pData_HighAcc$Distance < 700),]
  pData_LowAcc2 <- pData_LowAcc[(pData_LowAcc$Distance >= 350 & pData_LowAcc$Distance < 700),]
  pData_HighAcc3 <- pData_HighAcc[(pData_HighAcc$Distance >= 700 & pData_HighAcc$Distance < 1050),]
  pData_LowAcc3 <- pData_LowAcc[(pData_LowAcc$Distance >= 700 & pData_LowAcc$Distance < 1050),]
  pData_HighAcc4 <- pData_HighAcc[(pData_HighAcc$Distance >= 1050),]
  pData_LowAcc4 <- pData_LowAcc[(pData_LowAcc$Distance >= 1050),]
  
  # Statistics
  mean_pp[[p]] <- mean(pData_act2$ppLogNormalized)
  std_pp[[p]] <- sd(pData$ppLogNormalized)
  
  mean_pp_seg0[[p]] <- mean(pData_seg0$ppLogNormalized)
  mean_pp_seg0_1[[p]] <- mean(pData_seg0_1$ppLogNormalized)
  mean_pp_seg0_2[[p]] <- mean(pData_seg0_2$ppLogNormalized)
  mean_pp_seg0_3[[p]] <- mean(pData_seg0_3$ppLogNormalized)
  mean_pp_seg0_4[[p]] <- mean(pData_seg0_4$ppLogNormalized)
  
  mean_pp_seg1[[p]] <- mean(pData_seg1$ppLogNormalized)
  mean_pp_seg2[[p]] <- mean(pData_seg2$ppLogNormalized)
  mean_pp_seg3[[p]] <- mean(pData_seg3$ppLogNormalized)
  mean_pp_seg4[[p]] <- mean(pData_seg4$ppLogNormalized)
  mean_pp_max[[p]] <- max(mean_pp_seg1[[p]], mean_pp_seg2[[p]], mean_pp_seg3[[p]], mean_pp_seg4[[p]])
  
  std_pp_seg0[[p]] <- sd(pData_seg0$ppLogNormalized)
  std_pp_seg0_1[[p]] <- sd(pData_seg0_1$ppLogNormalized)
  std_pp_seg0_2[[p]] <- sd(pData_seg0_2$ppLogNormalized)
  std_pp_seg0_3[[p]] <- sd(pData_seg0_3$ppLogNormalized)
  std_pp_seg0_4[[p]] <- sd(pData_seg0_4$ppLogNormalized)
  
  std_pp_seg1[[p]] <- sd(pData_seg1$ppLogNormalized)
  std_pp_seg2[[p]] <- sd(pData_seg2$ppLogNormalized)
  std_pp_seg3[[p]] <- sd(pData_seg3$ppLogNormalized)
  std_pp_seg4[[p]] <- sd(pData_seg4$ppLogNormalized)
  std_pp_max[[p]] <- max(std_pp_seg1[[p]], std_pp_seg2[[p]], std_pp_seg3[[p]], std_pp_seg4[[p]])
  
  # Acceleration
  mean_pp_HighAcc[[p]] <- mean(pData_HighAcc$ppNext)
  mean_pp_LowAcc[[p]] <- mean(pData_LowAcc$ppNext)
  std_pp_HighAcc[[p]] <- sd(pData_HighAcc$ppNext)
  std_pp_LowAcc[[p]] <- sd(pData_LowAcc$ppNext)
  
  mean_pp_HighAcc1[[p]] <- mean(pData_HighAcc1$ppNext)
  mean_pp_LowAcc1[[p]] <- mean(pData_LowAcc1$ppNext)
  std_pp_HighAcc1[[p]] <- sd(pData_HighAcc1$ppNext)
  std_pp_LowAcc1[[p]] <- sd(pData_LowAcc1$ppNext)
  
  mean_pp_HighAcc2[[p]] <- mean(pData_HighAcc2$ppNext)
  mean_pp_LowAcc2[[p]] <- mean(pData_LowAcc2$ppNext)
  std_pp_HighAcc2[[p]] <- sd(pData_HighAcc2$ppNext)
  std_pp_LowAcc2[[p]] <- sd(pData_LowAcc2$ppNext)
  
  mean_pp_HighAcc3[[p]] <- mean(pData_HighAcc3$ppNext)
  mean_pp_LowAcc3[[p]] <- mean(pData_LowAcc3$ppNext)
  std_pp_HighAcc3[[p]] <- sd(pData_HighAcc3$ppNext)
  std_pp_LowAcc3[[p]] <- sd(pData_LowAcc3$ppNext)
  
  mean_pp_HighAcc4[[p]] <- mean(pData_HighAcc4$ppNext)
  mean_pp_LowAcc4[[p]] <- mean(pData_LowAcc4$ppNext)
  std_pp_HighAcc4[[p]] <- sd(pData_HighAcc4$ppNext)
  std_pp_LowAcc4[[p]] <- sd(pData_LowAcc4$ppNext)
}
plt_AllAcc <- vector(mode="list", length=length(persons)) 
names(plt_AllAcc) <- persons

COLOR_ACC = "#02A3C8"
COLOR_PP = "#F28E8E"
COLOR_BRAKE = "#888888"

y1 <- list(
  tickfont = list(color = COLOR_ACC),
  title="Degree",
  range=c(0, max(all_Drive2$Acceleration))
)
y2 <- list(
  tickfont = list(color = COLOR_PP),
  overlaying = "y",
  side = "right",
  title = "Log Perspiration",
  showgrid = FALSE,
  range=c(-0.6, 0.9)
  # range=c(min(all_Drive2$ppLogNormalized), max(all_Drive2$ppLogNormalized))
)

for (p in persons) {
  pData <- all_Drive2[(all_Drive2$Subject==as.integer(p) | all_Drive2$Subject==p),]
  
  pData_seg0 <- pData[pData$Phase==0,]
  pData_seg1 <- pData[pData$Phase==1 & pData$Activity==2 & pData$Time < 110,]
  pData_seg2 <- pData[pData$Phase==2 & pData$Activity==2 & pData$Time < 250,]
  pData_seg3 <- pData[pData$Phase==3 & pData$Activity==2 & pData$Time < 350,]
  pData_seg4 <- pData[pData$Phase==4 & pData$Activity==2,]
  
  plot_Acc <- plot_ly(pData, x = ~Time, height=400, width=900) %>%
    # add_trace(name="Acceleration", y = ~Acceleration, type = 'scatter', mode = 'lines', line=list(width=1.5, color=COLOR_ACC)) %>% 
    add_trace(name="PP", y = ~ppLogNormalized, type = 'scatter', mode = 'lines', connectgaps=F, line=list(width=1.5, color=COLOR_PP), yaxis = "y2") %>%
    add_segments(x = min(pData$Time), xend = max(pData$Time), y = mean_pp[[p]], yend = mean_pp[[p]],
                           yaxis = "y2", name="Avg. PP (straight)",
                           line=list(color="darkgray", dash = 'dot')) %>%
    add_segments(x = min(pData$Time), xend = max(pData$Time), y = mean_pp_seg0[[p]], yend = mean_pp_seg0[[p]], 
                           yaxis = "y2", name="Avg. PP (turning)",
                           line=list(color="black", dash = 'dot')) %>%
    add_segments(x = min(pData_seg1$Time), xend = max(pData_seg1$Time), y = mean_pp_seg1[[p]], yend = mean_pp_seg1[[p]], 
                           yaxis = "y2", name="Avg. PP (1st part)",
                           line=list(color="red", dash = 'dot')) %>%
    add_segments(x = min(pData_seg2$Time), xend = max(pData_seg2$Time), y = mean_pp_seg2[[p]], yend = mean_pp_seg2[[p]], 
                           yaxis = "y2", name="Avg. PP (2nd part)",
                           line=list(color="green", dash = 'dot')) %>%
    add_segments(x = min(pData_seg3$Time), xend = max(pData_seg3$Time), y = mean_pp_seg3[[p]], yend = mean_pp_seg3[[p]], 
                           yaxis = "y2", name="Avg. PP (3rd part)",
                           line=list(color="blue", dash = 'dot')) %>%
    add_segments(x = min(pData_seg4$Time), xend = max(pData_seg4$Time), y = mean_pp_seg4[[p]], yend = mean_pp_seg4[[p]], 
                           yaxis = "y2", name="Avg. PP (4th part)",
                           line=list(color="purple", dash = 'dot')) %>%
    layout(
      title=paste0("Subject #", p), 
      xaxis=list(title="Time [s]", range=c(0)), 
      yaxis=y1, 
      yaxis2=y2, 
      margin = list(l = 50, r = 50, b = 50, t = 50, pad = 4),
      legend = list(x = 0.5, xanchor = "center", y = 0.2, bgcolor = "rgba(0,0,0,0)", title="Metric", orientation = "h"),
      autosize = F
    )
  
  plt_AllAcc[[p]] <- plot_Acc
}


htmltools::tagList(plt_AllAcc)
NUMBER_OF_CLUSTERS = 3

color_darkpink = "#e75480"
CLUSTER_BRANCH_COLORS <- c("blue", "darkred", color_darkpink, "black")[1:NUMBER_OF_CLUSTERS]
CLUSTER_LABEL_COLORS <- c("blue", "darkred", color_darkpink, "black")[1:NUMBER_OF_CLUSTERS]


dfPP <- as.data.frame(cbind(
                            unlist(mean_pp), 
                            unlist(std_pp), 
                            unlist(mean_pp_seg0), 
                            unlist(mean_pp_seg1), 
                            unlist(mean_pp_seg2), 
                            unlist(mean_pp_seg3), 
                            unlist(mean_pp_seg4),
                            unlist(mean_pp_max),
                            unlist(std_pp_seg0), 
                            unlist(std_pp_seg1), 
                            unlist(std_pp_seg2), 
                            unlist(std_pp_seg3), 
                            unlist(std_pp_seg4),
                            unlist(std_pp_max),
                            unlist(mean_pp_seg0_1),
                            unlist(mean_pp_seg0_2),
                            unlist(mean_pp_seg0_3),
                            unlist(mean_pp_seg0_4),
                            unlist(std_pp_seg0_1),
                            unlist(std_pp_seg0_2),
                            unlist(std_pp_seg0_3),
                            unlist(std_pp_seg0_4),
                            unlist(mean_pp_HighAcc),
                            unlist(mean_pp_LowAcc),
                            unlist(std_pp_HighAcc),
                            unlist(std_pp_LowAcc),
                            unlist(mean_pp_HighAcc1),
                            unlist(mean_pp_LowAcc1),
                            unlist(std_pp_HighAcc1),
                            unlist(std_pp_LowAcc1),
                            unlist(mean_pp_HighAcc2),
                            unlist(mean_pp_LowAcc2),
                            unlist(std_pp_HighAcc2),
                            unlist(std_pp_LowAcc2),
                            unlist(mean_pp_HighAcc3),
                            unlist(mean_pp_LowAcc3),
                            unlist(std_pp_HighAcc3),
                            unlist(std_pp_LowAcc3),
                            unlist(mean_pp_HighAcc4),
                            unlist(mean_pp_LowAcc4),
                            unlist(std_pp_HighAcc4),
                            unlist(std_pp_LowAcc4)
              ))

names(dfPP) <- c("MeanPP", "StdPP", 
                 "MeanPP_Seg0", "MeanPP_Seg1", "MeanPP_Seg2", "MeanPP_Seg3", "MeanPP_Seg4", "MeanPP_SegMax",
                 "StdPP_Seg0", "StdPP_Seg1", "StdPP_Seg2", "StdPP_Seg3", "StdPP_Seg4", "StdPP_SegMax",
                 "MeanPP_Seg0_1", "MeanPP_Seg0_2", "MeanPP_Seg0_3", "MeanPP_Seg0_4",
                 "StdPP_Seg0_1", "StdPP_Seg0_2", "StdPP_Seg0_3", "StdPP_Seg0_4",
                 "MeanPP_AccHigh", " MeanPP_AccLow", "StdPP_AccHigh", "StdPP_AccLow",
                 "MeanPP_AccHigh1", " MeanPP_AccLow1", "StdPP_AccHigh1", "StdPP_AccLow1",
                 "MeanPP_AccHigh2", " MeanPP_AccLow2", "StdPP_AccHigh2", "StdPP_AccLow2",
                 "MeanPP_AccHigh3", " MeanPP_AccLow3", "StdPP_AccHigh3", "StdPP_AccLow3",
                 "MeanPP_AccHigh4", " MeanPP_AccLow4", "StdPP_AccHigh4", "StdPP_AccLow4")
behavioralMatrixClustering <- as.matrix(dfPP)

distMatrix <- dist(behavioralMatrixClustering)
hresults <- distMatrix %>% hclust

hc <- hresults %>% 
      as.dendrogram %>%
      set("nodes_cex", NUMBER_OF_CLUSTERS) %>%
      set("labels_col", value = CLUSTER_LABEL_COLORS, k=NUMBER_OF_CLUSTERS) %>%
      # set("leaves_pch", 19) %>%
      # set("leaves_col", value = c("gray"), k=NUMBER_OF_CLUSTERS) %>%    
      set("branches_k_color", value=CLUSTER_BRANCH_COLORS, k=NUMBER_OF_CLUSTERS)

plot(hc)
legend("topright", 
     title="Drive=Cognitive \nHierachical Clustering",
     legend = c("Group 1", "Group 2", "Group 3"),
     col = c("darkred", "pink" , "blue"),
     pch = c(20,20,20), bty = "n",  pt.cex = 1.5, cex = 0.8 , 
     text.col = "black", horiz = FALSE, inset = c(0.4, 0.1))

# Store clustering data
fPath <- str_interp("../data/processed/analysis/TT1_Drive_2_PP.csv")
dfx <- dfPP
dfx <- cbind(persons, dfx)
names(dfx) <- c("Subject", names(dfPP))
write.csv(dfx, fPath, row.names = F)
LS0tCnRpdGxlOiAiUiBOb3RlYm9vayIKb3V0cHV0OiBodG1sX25vdGVib29rCi0tLQoKYGBge3J9CnNvdXJjZSgnLi4vc2V0dGluZ3Mvc2V0dGluZ3MuUicpCnNvdXJjZSgnY29tbW9uRnVuY3Rpb25zLlInKQpgYGAKCmBgYHtyfQpwZXJzb25zIDwtIFNFTEVDVEVEX1NVQkpFQ1RTCgphbGxfRHJpdmUyIDwtIHJlYWQuY3N2KCcuLi8vZGF0YS9wcm9jZXNzZWQvZGlzdGFuY2V3aXNlL1RUMV9Ecml2ZV8yXzE1bV8xNW0uY3N2JykKYWxsX0RyaXZlMiRTdWJqZWN0IDwtIGFzLmZhY3RvcihhbGxfRHJpdmUyJFN1YmplY3QpCmFsbF9Ecml2ZTIkbG9nUGVyc3BpcmF0aW9uIDwtIGxvZyhhbGxfRHJpdmUyJFBlcnNwaXJhdGlvbikKYGBgCgojIEFjY2VsZXJhdGlvbiBQbG90CmBgYHtyfQpjb21wdXRlVGhyZXNob2xkIDwtIGZ1bmN0aW9uKHgsIHhtaW49bWluKHgpLCB4bWF4PW1heCh4KSkgewogIHggPC0geFt4ID49IHhtaW4gJiB4IDw9IHhtYXhdCiAgeEFyciA8LSBtYXRyaXgoeCwgbnJvdyA9IDEsIG5jb2wgPSBsZW5ndGgoeCkpCiAgdGhyZXNob2xkIDwtIG90c3UoeEFyciwgcmFuZ2U9Yyh4bWluLCB4bWF4KSkKICByZXR1cm4odGhyZXNob2xkKQp9CmNvbXB1dGVUaHJlc2hvbGRGb3JBY2NlbGVyYXRpb24gPC0gZnVuY3Rpb24oYWNjKXsKICByZXR1cm4oY29tcHV0ZVRocmVzaG9sZChhY2MsIHhtaW49MCwgeG1heD0yMCkpCn0KZmluZE1vZGVzIDwtIGZ1bmN0aW9uKHgpIHsKICB1eCA8LSB1bmlxdWUoeCkKICB1eFt3aGljaC5tYXgodGFidWxhdGUobWF0Y2goeCwgdXgpKSldCn0KYGBgCgpgYGB7cn0KYWNjX3RocmVzaG9sZHMgPC0gdmVjdG9yKG1vZGU9Imxpc3QiLCBsZW5ndGg9bGVuZ3RoKHBlcnNvbnMpKSAKbmFtZXMoYWNjX3RocmVzaG9sZHMpIDwtIHBlcnNvbnMKCmFjY19zZF90aHJlc2hvbGRzIDwtIHZlY3Rvcihtb2RlPSJsaXN0IiwgbGVuZ3RoPWxlbmd0aChwZXJzb25zKSkgCm5hbWVzKGFjY19zZF90aHJlc2hvbGRzKSA8LSBwZXJzb25zCgpzcGVlZF9zZF90aHJlc2hvbGRzIDwtIHZlY3Rvcihtb2RlPSJsaXN0IiwgbGVuZ3RoPWxlbmd0aChwZXJzb25zKSkgCm5hbWVzKHNwZWVkX3NkX3RocmVzaG9sZHMpIDwtIHBlcnNvbnMKCmZvciAocCBpbiBwZXJzb25zKSB7CiAgcERhdGEgPC0gYWxsX0RyaXZlMlsoYWxsX0RyaXZlMiRTdWJqZWN0PT1hcy5pbnRlZ2VyKHApIHwgYWxsX0RyaXZlMiRTdWJqZWN0PT1wKSxdCiAgYWNjX3RoIDwtIGNvbXB1dGVUaHJlc2hvbGRGb3JBY2NlbGVyYXRpb24ocERhdGEkQWNjZWxlcmF0aW9uKQogIGFjY190aF9zZCA8LSBjb21wdXRlVGhyZXNob2xkRm9yQWNjZWxlcmF0aW9uKHBEYXRhJEFjY19zdGQpCiAgc3BlZWRfdGhfc2QgPC0gY29tcHV0ZVRocmVzaG9sZEZvckFjY2VsZXJhdGlvbihwRGF0YSRTcGVlZF9zdGQpCiAgYWNjX3RocmVzaG9sZHNbW3BdXSA8LSBhcy5udW1lcmljKGFjY190aCkKICBhY2Nfc2RfdGhyZXNob2xkc1tbcF1dIDwtIGFzLm51bWVyaWMoYWNjX3RoX3NkKQogIHNwZWVkX3NkX3RocmVzaG9sZHNbW3BdXSA8LSBhcy5udW1lcmljKHNwZWVkX3RoX3NkKQogIAogIGFjY19wbCA8LSBnZ3Bsb3QocERhdGEsIGFlcyh4PUFjY2VsZXJhdGlvbikpICsgZ2VvbV9kZW5zaXR5KCkgKyBnZW9tX3ZsaW5lKGFlcyh4aW50ZXJjZXB0PWFjY190aCksIGNvbG9yPSJibHVlIiwgbGluZXR5cGU9ImRhc2hlZCIsIHNpemU9MSkKICBhY2Nfc2RfcGwgPC0gZ2dwbG90KHBEYXRhLCBhZXMoeD1BY2Nfc3RkKSkgKyBnZW9tX2RlbnNpdHkoKSArIGdlb21fdmxpbmUoYWVzKHhpbnRlcmNlcHQ9YWNjX3RoX3NkKSwgY29sb3I9InJlZCIsIGxpbmV0eXBlPSJkYXNoZWQiLCBzaXplPTEpCiAgc3BlZWRfc2RfcGwgPC0gZ2dwbG90KHBEYXRhLCBhZXMoeD1TcGVlZF9zdGQpKSArIGdlb21fZGVuc2l0eSgpICsgZ2VvbV92bGluZShhZXMoeGludGVyY2VwdD1zcGVlZF90aF9zZCksIGNvbG9yPSJyZWQiLCBsaW5ldHlwZT0iZGFzaGVkIiwgc2l6ZT0xKQogIAogIHBsb3QoYWNjX3BsKQogIHByaW50KHBhc3RlMChwLCAiIC0gVGhyZXNob2xkIChBY2MpID0iLCBhY2NfdGgsICIgJiAoQWNjIFNEKSA9IiwgYWNjX3RoX3NkLCAiICYgKFNwZWVkIFNEKSA9IiwgYWNjX3RoX3NkKSkKICAKfQpgYGAKCgpgYGB7cn0KbWVhbl9wcCA8LSB2ZWN0b3IobW9kZT0ibGlzdCIsIGxlbmd0aD1sZW5ndGgocGVyc29ucykpIApuYW1lcyhtZWFuX3BwKSA8LSBwZXJzb25zCgpzdGRfcHAgPC0gdmVjdG9yKG1vZGU9Imxpc3QiLCBsZW5ndGg9bGVuZ3RoKHBlcnNvbnMpKSAKbmFtZXMoc3RkX3BwKSA8LSBwZXJzb25zCgojIE1lYW4gKFR1cm5pbmcpCm1lYW5fcHBfc2VnMCA8LSB2ZWN0b3IobW9kZT0ibGlzdCIsIGxlbmd0aD1sZW5ndGgocGVyc29ucykpIApuYW1lcyhtZWFuX3BwX3NlZzApIDwtIHBlcnNvbnMKbWVhbl9wcF9zZWcwXzEgPC0gdmVjdG9yKG1vZGU9Imxpc3QiLCBsZW5ndGg9bGVuZ3RoKHBlcnNvbnMpKSAKbmFtZXMobWVhbl9wcF9zZWcwXzEpIDwtIHBlcnNvbnMKbWVhbl9wcF9zZWcwXzIgPC0gdmVjdG9yKG1vZGU9Imxpc3QiLCBsZW5ndGg9bGVuZ3RoKHBlcnNvbnMpKSAKbmFtZXMobWVhbl9wcF9zZWcwXzIpIDwtIHBlcnNvbnMKbWVhbl9wcF9zZWcwXzMgPC0gdmVjdG9yKG1vZGU9Imxpc3QiLCBsZW5ndGg9bGVuZ3RoKHBlcnNvbnMpKSAKbmFtZXMobWVhbl9wcF9zZWcwXzMpIDwtIHBlcnNvbnMKbWVhbl9wcF9zZWcwXzQgPC0gdmVjdG9yKG1vZGU9Imxpc3QiLCBsZW5ndGg9bGVuZ3RoKHBlcnNvbnMpKSAKbmFtZXMobWVhbl9wcF9zZWcwXzQpIDwtIHBlcnNvbnMKCiMgTWVhbiAoU3RyYWlnaHQpCm1lYW5fcHBfc2VnMSA8LSB2ZWN0b3IobW9kZT0ibGlzdCIsIGxlbmd0aD1sZW5ndGgocGVyc29ucykpIApuYW1lcyhtZWFuX3BwX3NlZzEpIDwtIHBlcnNvbnMKbWVhbl9wcF9zZWcyIDwtIHZlY3Rvcihtb2RlPSJsaXN0IiwgbGVuZ3RoPWxlbmd0aChwZXJzb25zKSkgCm5hbWVzKG1lYW5fcHBfc2VnMikgPC0gcGVyc29ucwptZWFuX3BwX3NlZzMgPC0gdmVjdG9yKG1vZGU9Imxpc3QiLCBsZW5ndGg9bGVuZ3RoKHBlcnNvbnMpKSAKbmFtZXMobWVhbl9wcF9zZWczKSA8LSBwZXJzb25zCm1lYW5fcHBfc2VnNCA8LSB2ZWN0b3IobW9kZT0ibGlzdCIsIGxlbmd0aD1sZW5ndGgocGVyc29ucykpIApuYW1lcyhtZWFuX3BwX3NlZzQpIDwtIHBlcnNvbnMKbWVhbl9wcF9tYXggPC0gdmVjdG9yKG1vZGU9Imxpc3QiLCBsZW5ndGg9bGVuZ3RoKHBlcnNvbnMpKSAKbmFtZXMobWVhbl9wcF9tYXgpIDwtIHBlcnNvbnMKCiMgU0QgKFR1cm5pbmcpCnN0ZF9wcF9zZWcwIDwtIHZlY3Rvcihtb2RlPSJsaXN0IiwgbGVuZ3RoPWxlbmd0aChwZXJzb25zKSkgCm5hbWVzKHN0ZF9wcF9zZWcwKSA8LSBwZXJzb25zCnN0ZF9wcF9zZWcwXzEgPC0gdmVjdG9yKG1vZGU9Imxpc3QiLCBsZW5ndGg9bGVuZ3RoKHBlcnNvbnMpKSAKbmFtZXMoc3RkX3BwX3NlZzBfMSkgPC0gcGVyc29ucwpzdGRfcHBfc2VnMF8yIDwtIHZlY3Rvcihtb2RlPSJsaXN0IiwgbGVuZ3RoPWxlbmd0aChwZXJzb25zKSkgCm5hbWVzKHN0ZF9wcF9zZWcwXzIpIDwtIHBlcnNvbnMKc3RkX3BwX3NlZzBfMyA8LSB2ZWN0b3IobW9kZT0ibGlzdCIsIGxlbmd0aD1sZW5ndGgocGVyc29ucykpIApuYW1lcyhzdGRfcHBfc2VnMF8zKSA8LSBwZXJzb25zCnN0ZF9wcF9zZWcwXzQgPC0gdmVjdG9yKG1vZGU9Imxpc3QiLCBsZW5ndGg9bGVuZ3RoKHBlcnNvbnMpKSAKbmFtZXMoc3RkX3BwX3NlZzBfNCkgPC0gcGVyc29ucwoKIyBTRCAoU3RyYWlnaHQpCnN0ZF9wcF9zZWcxIDwtIHZlY3Rvcihtb2RlPSJsaXN0IiwgbGVuZ3RoPWxlbmd0aChwZXJzb25zKSkgCm5hbWVzKHN0ZF9wcF9zZWcxKSA8LSBwZXJzb25zCnN0ZF9wcF9zZWcyIDwtIHZlY3Rvcihtb2RlPSJsaXN0IiwgbGVuZ3RoPWxlbmd0aChwZXJzb25zKSkgCm5hbWVzKHN0ZF9wcF9zZWcyKSA8LSBwZXJzb25zCnN0ZF9wcF9zZWczIDwtIHZlY3Rvcihtb2RlPSJsaXN0IiwgbGVuZ3RoPWxlbmd0aChwZXJzb25zKSkgCm5hbWVzKHN0ZF9wcF9zZWczKSA8LSBwZXJzb25zCnN0ZF9wcF9zZWc0IDwtIHZlY3Rvcihtb2RlPSJsaXN0IiwgbGVuZ3RoPWxlbmd0aChwZXJzb25zKSkgCm5hbWVzKHN0ZF9wcF9zZWc0KSA8LSBwZXJzb25zCnN0ZF9wcF9tYXggPC0gdmVjdG9yKG1vZGU9Imxpc3QiLCBsZW5ndGg9bGVuZ3RoKHBlcnNvbnMpKSAKbmFtZXMoc3RkX3BwX21heCkgPC0gcGVyc29ucwoKIyBBY2NlbGVyYXRpb24gUGFydGl0aW9uaW5nCm1lYW5fcHBfSGlnaEFjYyA8LSB2ZWN0b3IobW9kZT0ibGlzdCIsIGxlbmd0aD1sZW5ndGgocGVyc29ucykpIApuYW1lcyhtZWFuX3BwX0hpZ2hBY2MpIDwtIHBlcnNvbnMKbWVhbl9wcF9Mb3dBY2MgPC0gdmVjdG9yKG1vZGU9Imxpc3QiLCBsZW5ndGg9bGVuZ3RoKHBlcnNvbnMpKSAKbmFtZXMobWVhbl9wcF9Mb3dBY2MpIDwtIHBlcnNvbnMKCnN0ZF9wcF9IaWdoQWNjIDwtIHZlY3Rvcihtb2RlPSJsaXN0IiwgbGVuZ3RoPWxlbmd0aChwZXJzb25zKSkgCm5hbWVzKHN0ZF9wcF9IaWdoQWNjKSA8LSBwZXJzb25zCnN0ZF9wcF9Mb3dBY2MgPC0gdmVjdG9yKG1vZGU9Imxpc3QiLCBsZW5ndGg9bGVuZ3RoKHBlcnNvbnMpKSAKbmFtZXMoc3RkX3BwX0xvd0FjYykgPC0gcGVyc29ucwoKIyBBY2NlbCArIFNlZ21lbnRhdGlvbgptZWFuX3BwX0hpZ2hBY2MxIDwtIHZlY3Rvcihtb2RlPSJsaXN0IiwgbGVuZ3RoPWxlbmd0aChwZXJzb25zKSkgCm5hbWVzKG1lYW5fcHBfSGlnaEFjYzEpIDwtIHBlcnNvbnMKbWVhbl9wcF9Mb3dBY2MxIDwtIHZlY3Rvcihtb2RlPSJsaXN0IiwgbGVuZ3RoPWxlbmd0aChwZXJzb25zKSkgCm5hbWVzKG1lYW5fcHBfTG93QWNjMSkgPC0gcGVyc29ucwpzdGRfcHBfSGlnaEFjYzEgPC0gdmVjdG9yKG1vZGU9Imxpc3QiLCBsZW5ndGg9bGVuZ3RoKHBlcnNvbnMpKSAKbmFtZXMoc3RkX3BwX0hpZ2hBY2MxKSA8LSBwZXJzb25zCnN0ZF9wcF9Mb3dBY2MxIDwtIHZlY3Rvcihtb2RlPSJsaXN0IiwgbGVuZ3RoPWxlbmd0aChwZXJzb25zKSkgCm5hbWVzKHN0ZF9wcF9Mb3dBY2MxKSA8LSBwZXJzb25zCgptZWFuX3BwX0hpZ2hBY2MyIDwtIHZlY3Rvcihtb2RlPSJsaXN0IiwgbGVuZ3RoPWxlbmd0aChwZXJzb25zKSkgCm5hbWVzKG1lYW5fcHBfSGlnaEFjYzIpIDwtIHBlcnNvbnMKbWVhbl9wcF9Mb3dBY2MyIDwtIHZlY3Rvcihtb2RlPSJsaXN0IiwgbGVuZ3RoPWxlbmd0aChwZXJzb25zKSkgCm5hbWVzKG1lYW5fcHBfTG93QWNjMikgPC0gcGVyc29ucwpzdGRfcHBfSGlnaEFjYzIgPC0gdmVjdG9yKG1vZGU9Imxpc3QiLCBsZW5ndGg9bGVuZ3RoKHBlcnNvbnMpKSAKbmFtZXMoc3RkX3BwX0hpZ2hBY2MyKSA8LSBwZXJzb25zCnN0ZF9wcF9Mb3dBY2MyIDwtIHZlY3Rvcihtb2RlPSJsaXN0IiwgbGVuZ3RoPWxlbmd0aChwZXJzb25zKSkgCm5hbWVzKHN0ZF9wcF9Mb3dBY2MyKSA8LSBwZXJzb25zCgptZWFuX3BwX0hpZ2hBY2MzIDwtIHZlY3Rvcihtb2RlPSJsaXN0IiwgbGVuZ3RoPWxlbmd0aChwZXJzb25zKSkgCm5hbWVzKG1lYW5fcHBfSGlnaEFjYzMpIDwtIHBlcnNvbnMKbWVhbl9wcF9Mb3dBY2MzIDwtIHZlY3Rvcihtb2RlPSJsaXN0IiwgbGVuZ3RoPWxlbmd0aChwZXJzb25zKSkgCm5hbWVzKG1lYW5fcHBfTG93QWNjMykgPC0gcGVyc29ucwpzdGRfcHBfSGlnaEFjYzMgPC0gdmVjdG9yKG1vZGU9Imxpc3QiLCBsZW5ndGg9bGVuZ3RoKHBlcnNvbnMpKSAKbmFtZXMoc3RkX3BwX0hpZ2hBY2MzKSA8LSBwZXJzb25zCnN0ZF9wcF9Mb3dBY2MzIDwtIHZlY3Rvcihtb2RlPSJsaXN0IiwgbGVuZ3RoPWxlbmd0aChwZXJzb25zKSkgCm5hbWVzKHN0ZF9wcF9Mb3dBY2MzKSA8LSBwZXJzb25zCgptZWFuX3BwX0hpZ2hBY2M0IDwtIHZlY3Rvcihtb2RlPSJsaXN0IiwgbGVuZ3RoPWxlbmd0aChwZXJzb25zKSkgCm5hbWVzKG1lYW5fcHBfSGlnaEFjYzQpIDwtIHBlcnNvbnMKbWVhbl9wcF9Mb3dBY2M0IDwtIHZlY3Rvcihtb2RlPSJsaXN0IiwgbGVuZ3RoPWxlbmd0aChwZXJzb25zKSkgCm5hbWVzKG1lYW5fcHBfTG93QWNjNCkgPC0gcGVyc29ucwpzdGRfcHBfSGlnaEFjYzQgPC0gdmVjdG9yKG1vZGU9Imxpc3QiLCBsZW5ndGg9bGVuZ3RoKHBlcnNvbnMpKSAKbmFtZXMoc3RkX3BwX0hpZ2hBY2M0KSA8LSBwZXJzb25zCnN0ZF9wcF9Mb3dBY2M0IDwtIHZlY3Rvcihtb2RlPSJsaXN0IiwgbGVuZ3RoPWxlbmd0aChwZXJzb25zKSkgCm5hbWVzKHN0ZF9wcF9Mb3dBY2M0KSA8LSBwZXJzb25zCgoKZm9yKHAgaW4gcGVyc29ucykgewogIHBEYXRhIDwtIGFsbF9Ecml2ZTJbKGFsbF9Ecml2ZTIkU3ViamVjdD09YXMuaW50ZWdlcihwKSB8IGFsbF9Ecml2ZTIkU3ViamVjdD09cCksXQogIHBEYXRhX2FjdDIgPC0gcERhdGFbcERhdGEkQWN0aXZpdHk9PTIsXQogIAogICMgRGF0YSBQYXJ0aXRpb25pbmcKICBwRGF0YV9zZWcwIDwtIHBEYXRhW3BEYXRhJFBoYXNlPT0wLF0KICBwRGF0YV9zZWcwXzEgPC0gcERhdGFbcERhdGEkQWN0aXZpdHk9PTAgJiBwRGF0YSRQaGFzZT09MCAmIHBEYXRhJFRpbWUgPiAxMDAgJiBwRGF0YSRUaW1lIDwgMjAwLF0KICBwRGF0YV9zZWcwXzIgPC0gcERhdGFbcERhdGEkQWN0aXZpdHk9PTAgJiBwRGF0YSRQaGFzZT09MCAmIHBEYXRhJFRpbWUgPiAyMDAgJiBwRGF0YSRUaW1lIDwgMzAwLF0KICBwRGF0YV9zZWcwXzMgPC0gcERhdGFbcERhdGEkQWN0aXZpdHk9PTAgJiBwRGF0YSRQaGFzZT09MCAmIHBEYXRhJFRpbWUgPiAzMDAgJiBwRGF0YSRUaW1lIDwgNDAwLF0KICBwRGF0YV9zZWcwXzQgPC0gcERhdGFbcERhdGEkQWN0aXZpdHk9PTAgJiBwRGF0YSRQaGFzZT09MCAmIHBEYXRhJFRpbWUgPiA0MDAgJiBwRGF0YSRUaW1lIDwgNTAwLF0KICAKICBwRGF0YV9zZWcxIDwtIHBEYXRhW3BEYXRhJFBoYXNlPT0xICYgcERhdGEkQWN0aXZpdHk9PTIgJiBwRGF0YSRUaW1lIDwgMTEwLF0KICBwRGF0YV9zZWcyIDwtIHBEYXRhW3BEYXRhJFBoYXNlPT0yICYgcERhdGEkQWN0aXZpdHk9PTIgJiBwRGF0YSRUaW1lIDwgMjUwLF0KICBwRGF0YV9zZWczIDwtIHBEYXRhW3BEYXRhJFBoYXNlPT0zICYgcERhdGEkQWN0aXZpdHk9PTIgJiBwRGF0YSRUaW1lIDwgMzUwLF0KICBwRGF0YV9zZWc0IDwtIHBEYXRhW3BEYXRhJFBoYXNlPT00ICYgcERhdGEkQWN0aXZpdHk9PTIsXQogIAogIHBEYXRhX0hpZ2hBY2MgPC0gcERhdGFbcERhdGEkQWNjZWxlcmF0aW9uID49IGFjY190aHJlc2hvbGRzW1twXV0sXQogIHBEYXRhX0xvd0FjYyA8LSBwRGF0YVtwRGF0YSRBY2NlbGVyYXRpb24gPCBhY2NfdGhyZXNob2xkc1tbcF1dLF0KICAjIHBEYXRhX0hpZ2hBY2MgPC0gcERhdGFbcERhdGEkQWNjX3N0ZCA8IGFjY19zZF90aHJlc2hvbGRzW1twXV0gJiBwRGF0YSRTcGVlZF9zdGQgPCBzcGVlZF9zZF90aHJlc2hvbGRzW1twXV0sXQogICMgcERhdGFfTG93QWNjIDwtIHBEYXRhW3BEYXRhJEFjY19zdGQgPj0gYWNjX3NkX3RocmVzaG9sZHNbW3BdXSB8IHBEYXRhJFNwZWVkX3N0ZCA+PSBzcGVlZF9zZF90aHJlc2hvbGRzW1twXV0sXQogIAogIHBEYXRhX0hpZ2hBY2MxIDwtIHBEYXRhX0hpZ2hBY2NbKHBEYXRhX0hpZ2hBY2MkRGlzdGFuY2UgPCAzNTApLF0KICBwRGF0YV9Mb3dBY2MxIDwtIHBEYXRhX0xvd0FjY1socERhdGFfTG93QWNjJERpc3RhbmNlIDwgMzUwKSxdCiAgcERhdGFfSGlnaEFjYzIgPC0gcERhdGFfSGlnaEFjY1socERhdGFfSGlnaEFjYyREaXN0YW5jZSA+PSAzNTAgJiBwRGF0YV9IaWdoQWNjJERpc3RhbmNlIDwgNzAwKSxdCiAgcERhdGFfTG93QWNjMiA8LSBwRGF0YV9Mb3dBY2NbKHBEYXRhX0xvd0FjYyREaXN0YW5jZSA+PSAzNTAgJiBwRGF0YV9Mb3dBY2MkRGlzdGFuY2UgPCA3MDApLF0KICBwRGF0YV9IaWdoQWNjMyA8LSBwRGF0YV9IaWdoQWNjWyhwRGF0YV9IaWdoQWNjJERpc3RhbmNlID49IDcwMCAmIHBEYXRhX0hpZ2hBY2MkRGlzdGFuY2UgPCAxMDUwKSxdCiAgcERhdGFfTG93QWNjMyA8LSBwRGF0YV9Mb3dBY2NbKHBEYXRhX0xvd0FjYyREaXN0YW5jZSA+PSA3MDAgJiBwRGF0YV9Mb3dBY2MkRGlzdGFuY2UgPCAxMDUwKSxdCiAgcERhdGFfSGlnaEFjYzQgPC0gcERhdGFfSGlnaEFjY1socERhdGFfSGlnaEFjYyREaXN0YW5jZSA+PSAxMDUwKSxdCiAgcERhdGFfTG93QWNjNCA8LSBwRGF0YV9Mb3dBY2NbKHBEYXRhX0xvd0FjYyREaXN0YW5jZSA+PSAxMDUwKSxdCiAgCiAgIyBTdGF0aXN0aWNzCiAgbWVhbl9wcFtbcF1dIDwtIG1lYW4ocERhdGFfYWN0MiRwcExvZ05vcm1hbGl6ZWQpCiAgc3RkX3BwW1twXV0gPC0gc2QocERhdGEkcHBMb2dOb3JtYWxpemVkKQogIAogIG1lYW5fcHBfc2VnMFtbcF1dIDwtIG1lYW4ocERhdGFfc2VnMCRwcExvZ05vcm1hbGl6ZWQpCiAgbWVhbl9wcF9zZWcwXzFbW3BdXSA8LSBtZWFuKHBEYXRhX3NlZzBfMSRwcExvZ05vcm1hbGl6ZWQpCiAgbWVhbl9wcF9zZWcwXzJbW3BdXSA8LSBtZWFuKHBEYXRhX3NlZzBfMiRwcExvZ05vcm1hbGl6ZWQpCiAgbWVhbl9wcF9zZWcwXzNbW3BdXSA8LSBtZWFuKHBEYXRhX3NlZzBfMyRwcExvZ05vcm1hbGl6ZWQpCiAgbWVhbl9wcF9zZWcwXzRbW3BdXSA8LSBtZWFuKHBEYXRhX3NlZzBfNCRwcExvZ05vcm1hbGl6ZWQpCiAgCiAgbWVhbl9wcF9zZWcxW1twXV0gPC0gbWVhbihwRGF0YV9zZWcxJHBwTG9nTm9ybWFsaXplZCkKICBtZWFuX3BwX3NlZzJbW3BdXSA8LSBtZWFuKHBEYXRhX3NlZzIkcHBMb2dOb3JtYWxpemVkKQogIG1lYW5fcHBfc2VnM1tbcF1dIDwtIG1lYW4ocERhdGFfc2VnMyRwcExvZ05vcm1hbGl6ZWQpCiAgbWVhbl9wcF9zZWc0W1twXV0gPC0gbWVhbihwRGF0YV9zZWc0JHBwTG9nTm9ybWFsaXplZCkKICBtZWFuX3BwX21heFtbcF1dIDwtIG1heChtZWFuX3BwX3NlZzFbW3BdXSwgbWVhbl9wcF9zZWcyW1twXV0sIG1lYW5fcHBfc2VnM1tbcF1dLCBtZWFuX3BwX3NlZzRbW3BdXSkKICAKICBzdGRfcHBfc2VnMFtbcF1dIDwtIHNkKHBEYXRhX3NlZzAkcHBMb2dOb3JtYWxpemVkKQogIHN0ZF9wcF9zZWcwXzFbW3BdXSA8LSBzZChwRGF0YV9zZWcwXzEkcHBMb2dOb3JtYWxpemVkKQogIHN0ZF9wcF9zZWcwXzJbW3BdXSA8LSBzZChwRGF0YV9zZWcwXzIkcHBMb2dOb3JtYWxpemVkKQogIHN0ZF9wcF9zZWcwXzNbW3BdXSA8LSBzZChwRGF0YV9zZWcwXzMkcHBMb2dOb3JtYWxpemVkKQogIHN0ZF9wcF9zZWcwXzRbW3BdXSA8LSBzZChwRGF0YV9zZWcwXzQkcHBMb2dOb3JtYWxpemVkKQogIAogIHN0ZF9wcF9zZWcxW1twXV0gPC0gc2QocERhdGFfc2VnMSRwcExvZ05vcm1hbGl6ZWQpCiAgc3RkX3BwX3NlZzJbW3BdXSA8LSBzZChwRGF0YV9zZWcyJHBwTG9nTm9ybWFsaXplZCkKICBzdGRfcHBfc2VnM1tbcF1dIDwtIHNkKHBEYXRhX3NlZzMkcHBMb2dOb3JtYWxpemVkKQogIHN0ZF9wcF9zZWc0W1twXV0gPC0gc2QocERhdGFfc2VnNCRwcExvZ05vcm1hbGl6ZWQpCiAgc3RkX3BwX21heFtbcF1dIDwtIG1heChzdGRfcHBfc2VnMVtbcF1dLCBzdGRfcHBfc2VnMltbcF1dLCBzdGRfcHBfc2VnM1tbcF1dLCBzdGRfcHBfc2VnNFtbcF1dKQogIAogICMgQWNjZWxlcmF0aW9uCiAgbWVhbl9wcF9IaWdoQWNjW1twXV0gPC0gbWVhbihwRGF0YV9IaWdoQWNjJHBwTmV4dCkKICBtZWFuX3BwX0xvd0FjY1tbcF1dIDwtIG1lYW4ocERhdGFfTG93QWNjJHBwTmV4dCkKICBzdGRfcHBfSGlnaEFjY1tbcF1dIDwtIHNkKHBEYXRhX0hpZ2hBY2MkcHBOZXh0KQogIHN0ZF9wcF9Mb3dBY2NbW3BdXSA8LSBzZChwRGF0YV9Mb3dBY2MkcHBOZXh0KQogIAogIG1lYW5fcHBfSGlnaEFjYzFbW3BdXSA8LSBtZWFuKHBEYXRhX0hpZ2hBY2MxJHBwTmV4dCkKICBtZWFuX3BwX0xvd0FjYzFbW3BdXSA8LSBtZWFuKHBEYXRhX0xvd0FjYzEkcHBOZXh0KQogIHN0ZF9wcF9IaWdoQWNjMVtbcF1dIDwtIHNkKHBEYXRhX0hpZ2hBY2MxJHBwTmV4dCkKICBzdGRfcHBfTG93QWNjMVtbcF1dIDwtIHNkKHBEYXRhX0xvd0FjYzEkcHBOZXh0KQogIAogIG1lYW5fcHBfSGlnaEFjYzJbW3BdXSA8LSBtZWFuKHBEYXRhX0hpZ2hBY2MyJHBwTmV4dCkKICBtZWFuX3BwX0xvd0FjYzJbW3BdXSA8LSBtZWFuKHBEYXRhX0xvd0FjYzIkcHBOZXh0KQogIHN0ZF9wcF9IaWdoQWNjMltbcF1dIDwtIHNkKHBEYXRhX0hpZ2hBY2MyJHBwTmV4dCkKICBzdGRfcHBfTG93QWNjMltbcF1dIDwtIHNkKHBEYXRhX0xvd0FjYzIkcHBOZXh0KQogIAogIG1lYW5fcHBfSGlnaEFjYzNbW3BdXSA8LSBtZWFuKHBEYXRhX0hpZ2hBY2MzJHBwTmV4dCkKICBtZWFuX3BwX0xvd0FjYzNbW3BdXSA8LSBtZWFuKHBEYXRhX0xvd0FjYzMkcHBOZXh0KQogIHN0ZF9wcF9IaWdoQWNjM1tbcF1dIDwtIHNkKHBEYXRhX0hpZ2hBY2MzJHBwTmV4dCkKICBzdGRfcHBfTG93QWNjM1tbcF1dIDwtIHNkKHBEYXRhX0xvd0FjYzMkcHBOZXh0KQogIAogIG1lYW5fcHBfSGlnaEFjYzRbW3BdXSA8LSBtZWFuKHBEYXRhX0hpZ2hBY2M0JHBwTmV4dCkKICBtZWFuX3BwX0xvd0FjYzRbW3BdXSA8LSBtZWFuKHBEYXRhX0xvd0FjYzQkcHBOZXh0KQogIHN0ZF9wcF9IaWdoQWNjNFtbcF1dIDwtIHNkKHBEYXRhX0hpZ2hBY2M0JHBwTmV4dCkKICBzdGRfcHBfTG93QWNjNFtbcF1dIDwtIHNkKHBEYXRhX0xvd0FjYzQkcHBOZXh0KQp9CgpgYGAKCmBgYHtyfQpwbHRfQWxsQWNjIDwtIHZlY3Rvcihtb2RlPSJsaXN0IiwgbGVuZ3RoPWxlbmd0aChwZXJzb25zKSkgCm5hbWVzKHBsdF9BbGxBY2MpIDwtIHBlcnNvbnMKCkNPTE9SX0FDQyA9ICIjMDJBM0M4IgpDT0xPUl9QUCA9ICIjRjI4RThFIgpDT0xPUl9CUkFLRSA9ICIjODg4ODg4IgoKeTEgPC0gbGlzdCgKICB0aWNrZm9udCA9IGxpc3QoY29sb3IgPSBDT0xPUl9BQ0MpLAogIHRpdGxlPSJEZWdyZWUiLAogIHJhbmdlPWMoMCwgbWF4KGFsbF9Ecml2ZTIkQWNjZWxlcmF0aW9uKSkKKQp5MiA8LSBsaXN0KAogIHRpY2tmb250ID0gbGlzdChjb2xvciA9IENPTE9SX1BQKSwKICBvdmVybGF5aW5nID0gInkiLAogIHNpZGUgPSAicmlnaHQiLAogIHRpdGxlID0gIkxvZyBQZXJzcGlyYXRpb24iLAogIHNob3dncmlkID0gRkFMU0UsCiAgcmFuZ2U9YygtMC42LCAwLjkpCiAgIyByYW5nZT1jKG1pbihhbGxfRHJpdmUyJHBwTG9nTm9ybWFsaXplZCksIG1heChhbGxfRHJpdmUyJHBwTG9nTm9ybWFsaXplZCkpCikKCmZvciAocCBpbiBwZXJzb25zKSB7CiAgcERhdGEgPC0gYWxsX0RyaXZlMlsoYWxsX0RyaXZlMiRTdWJqZWN0PT1hcy5pbnRlZ2VyKHApIHwgYWxsX0RyaXZlMiRTdWJqZWN0PT1wKSxdCiAgCiAgcERhdGFfc2VnMCA8LSBwRGF0YVtwRGF0YSRQaGFzZT09MCxdCiAgcERhdGFfc2VnMSA8LSBwRGF0YVtwRGF0YSRQaGFzZT09MSAmIHBEYXRhJEFjdGl2aXR5PT0yICYgcERhdGEkVGltZSA8IDExMCxdCiAgcERhdGFfc2VnMiA8LSBwRGF0YVtwRGF0YSRQaGFzZT09MiAmIHBEYXRhJEFjdGl2aXR5PT0yICYgcERhdGEkVGltZSA8IDI1MCxdCiAgcERhdGFfc2VnMyA8LSBwRGF0YVtwRGF0YSRQaGFzZT09MyAmIHBEYXRhJEFjdGl2aXR5PT0yICYgcERhdGEkVGltZSA8IDM1MCxdCiAgcERhdGFfc2VnNCA8LSBwRGF0YVtwRGF0YSRQaGFzZT09NCAmIHBEYXRhJEFjdGl2aXR5PT0yLF0KICAKICBwbG90X0FjYyA8LSBwbG90X2x5KHBEYXRhLCB4ID0gflRpbWUsIGhlaWdodD00MDAsIHdpZHRoPTkwMCkgJT4lCiAgICAjIGFkZF90cmFjZShuYW1lPSJBY2NlbGVyYXRpb24iLCB5ID0gfkFjY2VsZXJhdGlvbiwgdHlwZSA9ICdzY2F0dGVyJywgbW9kZSA9ICdsaW5lcycsIGxpbmU9bGlzdCh3aWR0aD0xLjUsIGNvbG9yPUNPTE9SX0FDQykpICU+JSAKICAgIGFkZF90cmFjZShuYW1lPSJQUCIsIHkgPSB+cHBMb2dOb3JtYWxpemVkLCB0eXBlID0gJ3NjYXR0ZXInLCBtb2RlID0gJ2xpbmVzJywgY29ubmVjdGdhcHM9RiwgbGluZT1saXN0KHdpZHRoPTEuNSwgY29sb3I9Q09MT1JfUFApLCB5YXhpcyA9ICJ5MiIpICU+JQogICAgYWRkX3NlZ21lbnRzKHggPSBtaW4ocERhdGEkVGltZSksIHhlbmQgPSBtYXgocERhdGEkVGltZSksIHkgPSBtZWFuX3BwW1twXV0sIHllbmQgPSBtZWFuX3BwW1twXV0sCiAgICAgICAgICAgICAgICAgICAgICAgICAgIHlheGlzID0gInkyIiwgbmFtZT0iQXZnLiBQUCAoc3RyYWlnaHQpIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgbGluZT1saXN0KGNvbG9yPSJkYXJrZ3JheSIsIGRhc2ggPSAnZG90JykpICU+JQogICAgYWRkX3NlZ21lbnRzKHggPSBtaW4ocERhdGEkVGltZSksIHhlbmQgPSBtYXgocERhdGEkVGltZSksIHkgPSBtZWFuX3BwX3NlZzBbW3BdXSwgeWVuZCA9IG1lYW5fcHBfc2VnMFtbcF1dLCAKICAgICAgICAgICAgICAgICAgICAgICAgICAgeWF4aXMgPSAieTIiLCBuYW1lPSJBdmcuIFBQICh0dXJuaW5nKSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgIGxpbmU9bGlzdChjb2xvcj0iYmxhY2siLCBkYXNoID0gJ2RvdCcpKSAlPiUKICAgIGFkZF9zZWdtZW50cyh4ID0gbWluKHBEYXRhX3NlZzEkVGltZSksIHhlbmQgPSBtYXgocERhdGFfc2VnMSRUaW1lKSwgeSA9IG1lYW5fcHBfc2VnMVtbcF1dLCB5ZW5kID0gbWVhbl9wcF9zZWcxW1twXV0sIAogICAgICAgICAgICAgICAgICAgICAgICAgICB5YXhpcyA9ICJ5MiIsIG5hbWU9IkF2Zy4gUFAgKDFzdCBwYXJ0KSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgIGxpbmU9bGlzdChjb2xvcj0icmVkIiwgZGFzaCA9ICdkb3QnKSkgJT4lCiAgICBhZGRfc2VnbWVudHMoeCA9IG1pbihwRGF0YV9zZWcyJFRpbWUpLCB4ZW5kID0gbWF4KHBEYXRhX3NlZzIkVGltZSksIHkgPSBtZWFuX3BwX3NlZzJbW3BdXSwgeWVuZCA9IG1lYW5fcHBfc2VnMltbcF1dLCAKICAgICAgICAgICAgICAgICAgICAgICAgICAgeWF4aXMgPSAieTIiLCBuYW1lPSJBdmcuIFBQICgybmQgcGFydCkiLAogICAgICAgICAgICAgICAgICAgICAgICAgICBsaW5lPWxpc3QoY29sb3I9ImdyZWVuIiwgZGFzaCA9ICdkb3QnKSkgJT4lCiAgICBhZGRfc2VnbWVudHMoeCA9IG1pbihwRGF0YV9zZWczJFRpbWUpLCB4ZW5kID0gbWF4KHBEYXRhX3NlZzMkVGltZSksIHkgPSBtZWFuX3BwX3NlZzNbW3BdXSwgeWVuZCA9IG1lYW5fcHBfc2VnM1tbcF1dLCAKICAgICAgICAgICAgICAgICAgICAgICAgICAgeWF4aXMgPSAieTIiLCBuYW1lPSJBdmcuIFBQICgzcmQgcGFydCkiLAogICAgICAgICAgICAgICAgICAgICAgICAgICBsaW5lPWxpc3QoY29sb3I9ImJsdWUiLCBkYXNoID0gJ2RvdCcpKSAlPiUKICAgIGFkZF9zZWdtZW50cyh4ID0gbWluKHBEYXRhX3NlZzQkVGltZSksIHhlbmQgPSBtYXgocERhdGFfc2VnNCRUaW1lKSwgeSA9IG1lYW5fcHBfc2VnNFtbcF1dLCB5ZW5kID0gbWVhbl9wcF9zZWc0W1twXV0sIAogICAgICAgICAgICAgICAgICAgICAgICAgICB5YXhpcyA9ICJ5MiIsIG5hbWU9IkF2Zy4gUFAgKDR0aCBwYXJ0KSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgIGxpbmU9bGlzdChjb2xvcj0icHVycGxlIiwgZGFzaCA9ICdkb3QnKSkgJT4lCiAgICBsYXlvdXQoCiAgICAgIHRpdGxlPXBhc3RlMCgiU3ViamVjdCAjIiwgcCksIAogICAgICB4YXhpcz1saXN0KHRpdGxlPSJUaW1lIFtzXSIsIHJhbmdlPWMoMCkpLCAKICAgICAgeWF4aXM9eTEsIAogICAgICB5YXhpczI9eTIsIAogICAgICBtYXJnaW4gPSBsaXN0KGwgPSA1MCwgciA9IDUwLCBiID0gNTAsIHQgPSA1MCwgcGFkID0gNCksCiAgICAgIGxlZ2VuZCA9IGxpc3QoeCA9IDAuNSwgeGFuY2hvciA9ICJjZW50ZXIiLCB5ID0gMC4yLCBiZ2NvbG9yID0gInJnYmEoMCwwLDAsMCkiLCB0aXRsZT0iTWV0cmljIiwgb3JpZW50YXRpb24gPSAiaCIpLAogICAgICBhdXRvc2l6ZSA9IEYKICAgICkKICAKICBwbHRfQWxsQWNjW1twXV0gPC0gcGxvdF9BY2MKfQoKCmh0bWx0b29sczo6dGFnTGlzdChwbHRfQWxsQWNjKQpgYGAKCgpgYGB7cn0KTlVNQkVSX09GX0NMVVNURVJTID0gMwoKY29sb3JfZGFya3BpbmsgPSAiI2U3NTQ4MCIKQ0xVU1RFUl9CUkFOQ0hfQ09MT1JTIDwtIGMoImJsdWUiLCAiZGFya3JlZCIsIGNvbG9yX2RhcmtwaW5rLCAiYmxhY2siKVsxOk5VTUJFUl9PRl9DTFVTVEVSU10KQ0xVU1RFUl9MQUJFTF9DT0xPUlMgPC0gYygiYmx1ZSIsICJkYXJrcmVkIiwgY29sb3JfZGFya3BpbmssICJibGFjayIpWzE6TlVNQkVSX09GX0NMVVNURVJTXQoKCmRmUFAgPC0gYXMuZGF0YS5mcmFtZShjYmluZCgKICAgICAgICAgICAgICAgICAgICAgICAgICAgIHVubGlzdChtZWFuX3BwKSwgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICB1bmxpc3Qoc3RkX3BwKSwgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICB1bmxpc3QobWVhbl9wcF9zZWcwKSwgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICB1bmxpc3QobWVhbl9wcF9zZWcxKSwgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICB1bmxpc3QobWVhbl9wcF9zZWcyKSwgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICB1bmxpc3QobWVhbl9wcF9zZWczKSwgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICB1bmxpc3QobWVhbl9wcF9zZWc0KSwKICAgICAgICAgICAgICAgICAgICAgICAgICAgIHVubGlzdChtZWFuX3BwX21heCksCiAgICAgICAgICAgICAgICAgICAgICAgICAgICB1bmxpc3Qoc3RkX3BwX3NlZzApLCAKICAgICAgICAgICAgICAgICAgICAgICAgICAgIHVubGlzdChzdGRfcHBfc2VnMSksIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgdW5saXN0KHN0ZF9wcF9zZWcyKSwgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICB1bmxpc3Qoc3RkX3BwX3NlZzMpLCAKICAgICAgICAgICAgICAgICAgICAgICAgICAgIHVubGlzdChzdGRfcHBfc2VnNCksCiAgICAgICAgICAgICAgICAgICAgICAgICAgICB1bmxpc3Qoc3RkX3BwX21heCksCiAgICAgICAgICAgICAgICAgICAgICAgICAgICB1bmxpc3QobWVhbl9wcF9zZWcwXzEpLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgdW5saXN0KG1lYW5fcHBfc2VnMF8yKSwKICAgICAgICAgICAgICAgICAgICAgICAgICAgIHVubGlzdChtZWFuX3BwX3NlZzBfMyksCiAgICAgICAgICAgICAgICAgICAgICAgICAgICB1bmxpc3QobWVhbl9wcF9zZWcwXzQpLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgdW5saXN0KHN0ZF9wcF9zZWcwXzEpLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgdW5saXN0KHN0ZF9wcF9zZWcwXzIpLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgdW5saXN0KHN0ZF9wcF9zZWcwXzMpLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgdW5saXN0KHN0ZF9wcF9zZWcwXzQpLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgdW5saXN0KG1lYW5fcHBfSGlnaEFjYyksCiAgICAgICAgICAgICAgICAgICAgICAgICAgICB1bmxpc3QobWVhbl9wcF9Mb3dBY2MpLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgdW5saXN0KHN0ZF9wcF9IaWdoQWNjKSwKICAgICAgICAgICAgICAgICAgICAgICAgICAgIHVubGlzdChzdGRfcHBfTG93QWNjKSwKICAgICAgICAgICAgICAgICAgICAgICAgICAgIHVubGlzdChtZWFuX3BwX0hpZ2hBY2MxKSwKICAgICAgICAgICAgICAgICAgICAgICAgICAgIHVubGlzdChtZWFuX3BwX0xvd0FjYzEpLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgdW5saXN0KHN0ZF9wcF9IaWdoQWNjMSksCiAgICAgICAgICAgICAgICAgICAgICAgICAgICB1bmxpc3Qoc3RkX3BwX0xvd0FjYzEpLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgdW5saXN0KG1lYW5fcHBfSGlnaEFjYzIpLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgdW5saXN0KG1lYW5fcHBfTG93QWNjMiksCiAgICAgICAgICAgICAgICAgICAgICAgICAgICB1bmxpc3Qoc3RkX3BwX0hpZ2hBY2MyKSwKICAgICAgICAgICAgICAgICAgICAgICAgICAgIHVubGlzdChzdGRfcHBfTG93QWNjMiksCiAgICAgICAgICAgICAgICAgICAgICAgICAgICB1bmxpc3QobWVhbl9wcF9IaWdoQWNjMyksCiAgICAgICAgICAgICAgICAgICAgICAgICAgICB1bmxpc3QobWVhbl9wcF9Mb3dBY2MzKSwKICAgICAgICAgICAgICAgICAgICAgICAgICAgIHVubGlzdChzdGRfcHBfSGlnaEFjYzMpLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgdW5saXN0KHN0ZF9wcF9Mb3dBY2MzKSwKICAgICAgICAgICAgICAgICAgICAgICAgICAgIHVubGlzdChtZWFuX3BwX0hpZ2hBY2M0KSwKICAgICAgICAgICAgICAgICAgICAgICAgICAgIHVubGlzdChtZWFuX3BwX0xvd0FjYzQpLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgdW5saXN0KHN0ZF9wcF9IaWdoQWNjNCksCiAgICAgICAgICAgICAgICAgICAgICAgICAgICB1bmxpc3Qoc3RkX3BwX0xvd0FjYzQpCiAgICAgICAgICAgICAgKSkKCm5hbWVzKGRmUFApIDwtIGMoIk1lYW5QUCIsICJTdGRQUCIsIAogICAgICAgICAgICAgICAgICJNZWFuUFBfU2VnMCIsICJNZWFuUFBfU2VnMSIsICJNZWFuUFBfU2VnMiIsICJNZWFuUFBfU2VnMyIsICJNZWFuUFBfU2VnNCIsICJNZWFuUFBfU2VnTWF4IiwKICAgICAgICAgICAgICAgICAiU3RkUFBfU2VnMCIsICJTdGRQUF9TZWcxIiwgIlN0ZFBQX1NlZzIiLCAiU3RkUFBfU2VnMyIsICJTdGRQUF9TZWc0IiwgIlN0ZFBQX1NlZ01heCIsCiAgICAgICAgICAgICAgICAgIk1lYW5QUF9TZWcwXzEiLCAiTWVhblBQX1NlZzBfMiIsICJNZWFuUFBfU2VnMF8zIiwgIk1lYW5QUF9TZWcwXzQiLAogICAgICAgICAgICAgICAgICJTdGRQUF9TZWcwXzEiLCAiU3RkUFBfU2VnMF8yIiwgIlN0ZFBQX1NlZzBfMyIsICJTdGRQUF9TZWcwXzQiLAogICAgICAgICAgICAgICAgICJNZWFuUFBfQWNjSGlnaCIsICIgTWVhblBQX0FjY0xvdyIsICJTdGRQUF9BY2NIaWdoIiwgIlN0ZFBQX0FjY0xvdyIsCiAgICAgICAgICAgICAgICAgIk1lYW5QUF9BY2NIaWdoMSIsICIgTWVhblBQX0FjY0xvdzEiLCAiU3RkUFBfQWNjSGlnaDEiLCAiU3RkUFBfQWNjTG93MSIsCiAgICAgICAgICAgICAgICAgIk1lYW5QUF9BY2NIaWdoMiIsICIgTWVhblBQX0FjY0xvdzIiLCAiU3RkUFBfQWNjSGlnaDIiLCAiU3RkUFBfQWNjTG93MiIsCiAgICAgICAgICAgICAgICAgIk1lYW5QUF9BY2NIaWdoMyIsICIgTWVhblBQX0FjY0xvdzMiLCAiU3RkUFBfQWNjSGlnaDMiLCAiU3RkUFBfQWNjTG93MyIsCiAgICAgICAgICAgICAgICAgIk1lYW5QUF9BY2NIaWdoNCIsICIgTWVhblBQX0FjY0xvdzQiLCAiU3RkUFBfQWNjSGlnaDQiLCAiU3RkUFBfQWNjTG93NCIpCmJlaGF2aW9yYWxNYXRyaXhDbHVzdGVyaW5nIDwtIGFzLm1hdHJpeChkZlBQKQoKZGlzdE1hdHJpeCA8LSBkaXN0KGJlaGF2aW9yYWxNYXRyaXhDbHVzdGVyaW5nKQpocmVzdWx0cyA8LSBkaXN0TWF0cml4ICU+JSBoY2x1c3QKCmhjIDwtIGhyZXN1bHRzICU+JSAKICAgICAgYXMuZGVuZHJvZ3JhbSAlPiUKICAgICAgc2V0KCJub2Rlc19jZXgiLCBOVU1CRVJfT0ZfQ0xVU1RFUlMpICU+JQogICAgICBzZXQoImxhYmVsc19jb2wiLCB2YWx1ZSA9IENMVVNURVJfTEFCRUxfQ09MT1JTLCBrPU5VTUJFUl9PRl9DTFVTVEVSUykgJT4lCiAgICAgICMgc2V0KCJsZWF2ZXNfcGNoIiwgMTkpICU+JQogICAgICAjIHNldCgibGVhdmVzX2NvbCIsIHZhbHVlID0gYygiZ3JheSIpLCBrPU5VTUJFUl9PRl9DTFVTVEVSUykgJT4lICAgIAogICAgICBzZXQoImJyYW5jaGVzX2tfY29sb3IiLCB2YWx1ZT1DTFVTVEVSX0JSQU5DSF9DT0xPUlMsIGs9TlVNQkVSX09GX0NMVVNURVJTKQoKcGxvdChoYykKbGVnZW5kKCJ0b3ByaWdodCIsIAogICAgIHRpdGxlPSJEcml2ZT1Db2duaXRpdmUgXG5IaWVyYWNoaWNhbCBDbHVzdGVyaW5nIiwKICAgICBsZWdlbmQgPSBjKCJHcm91cCAxIiwgIkdyb3VwIDIiLCAiR3JvdXAgMyIpLAogICAgIGNvbCA9IGMoImRhcmtyZWQiLCAicGluayIgLCAiYmx1ZSIpLAogICAgIHBjaCA9IGMoMjAsMjAsMjApLCBidHkgPSAibiIsICBwdC5jZXggPSAxLjUsIGNleCA9IDAuOCAsIAogICAgIHRleHQuY29sID0gImJsYWNrIiwgaG9yaXogPSBGQUxTRSwgaW5zZXQgPSBjKDAuNCwgMC4xKSkKYGBgCgoKCmBgYHtyfQojIFN0b3JlIGNsdXN0ZXJpbmcgZGF0YQpmUGF0aCA8LSBzdHJfaW50ZXJwKCIuLi9kYXRhL3Byb2Nlc3NlZC9hbmFseXNpcy9UVDFfRHJpdmVfMl9QUC5jc3YiKQpkZnggPC0gZGZQUApkZnggPC0gY2JpbmQocGVyc29ucywgZGZ4KQpuYW1lcyhkZngpIDwtIGMoIlN1YmplY3QiLCBuYW1lcyhkZlBQKSkKd3JpdGUuY3N2KGRmeCwgZlBhdGgsIHJvdy5uYW1lcyA9IEYpCmBgYAoKCgo=